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ABSTRACT:  The  effects  of  the  injection  of  a  light,  medium,  and 
heavy  gas  on  the  skin  friction  and  heat  transfer  characteristics 
of  the  laminar  boundary  layer  flow  over  a  flat  plate  have  been 
investigated.  This  was  done  by  a  numerical  solution  of  the 
governing  differential  equations  for  injection  rate  parameters  of 
0,  -0.1,  -0.2;  Mach  numbers  from  0  to  3;  and  wall-to-stream  tem¬ 
perature  ratios  of  0.25,  0.50,  and  1.00.  The  results  show  that 
the  selection  of  discrete  injection  systems  must  represent  a 
favorable  compromise  between  conflicting  factors  (such  as  small 
molecular  weight  and  large  molecular  diameters)  which  influence 
the  amount  of  reduction  in  skin  friction  and  heat  transfer. 
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Skin  Friction  and  Heat  Transfer  Characteristics  of  the  Compressible 
Laminar  Boundary  Layer  with  Injection  of  a  Light,  Medium,  and  Heavy 
Gas 

This  is  the  first  part  of  a  three-fold  investigation  into  the 
effects  of  foreign  gas  injection  on  the  laminar  boundary  layer 
heat  transfer  and  stability  characteristics.  In  this  report  the 
boundary  layer  equations  are  solved  numerically  and  profiles  are 
obtained  showing  the  effects  of  light,  medium,  and  heavy  gas 
injection.  In  addition,  the  effects  of  the  injection  on  the  skin 
friction  and  heat  transfer  characteristics  of  these  profiles  are 
examined. 
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Reynolds  number  based  on  main  stream  conditions 
stagnation  enthalpy  (c.f.,  Eq.  (25)) 
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SUBSCRIPTS 

1  refers  to  injected  gas 

2  refers  to  main  stream  gas 
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INTRODUCTION 

This  report  describes  the  skin  friction  and  heat  transfer 
investigations  which  were  part  of  a  study  of  the  effects  of  foreign 
gas  injection  on  the  heat  transfer  and  stability  characteristics  of 
the  lasinar  boundary  layer  flow  over  a  flat  plate.  Specifically, 
the  effect  of  molecular  properties  of  the  injected  gas  on  the  skin 
friction  and  heat  transfer  is  considered  here.  The  stability 
characteristics  of  sose  of  the  resulting  boundary  layer  profiles 
appear  in  reference  (1) .  In  addition,  an  outline  of  a  numerical 
method  used  for  a  direct  numerical  solution  of  the  complete  system 
of  disturbance  equations  for  the  boundary  layer  stability  analysis 
is  given  in  reference  (2) . 

The  effects  of  injection  on  the  skin  friction  and  heat  transfer 
are  of  Interest  in  the  area  of  high  speed  vehicle  design,  where 
injection  of  a  foreign  gas  through  the  porous  wall  of  the  body  is 
recognized  as  an  effective  way  of  alleviating  surface  heat  transfer. 
The  present  investigation  is  based  on  the  analysis  of  Sken  (ref.  (3)> 
w!.o  derived  a  similar  ity  form  of  the  mass  transfer  boundary  layer 
equations  by  assuming  that  the  injection  velocity  varied  inversely 
with  the  square  root  of  the  distance  from  the  leading  edge  of  a  flat 
plate.  In  his  report,  Shen  did  not  include  any  numerical  results. 

However,  Korobkin,  in  reference  (4),  did  give  some  results 
showing  that  the  laminar  boundary  layer  heat  transfer  characteristics 
can  be  appreciably  altered  depending  on  the  molecular  properties  of 
the  injected  gas.  Specifically,  his  analysis  demonstrated  that 
small  molecular  weight  and  large  molecular  diameter  of  the  injected 
gas  are  desirable  properties  for  the  reduction  of  skin  friction  and 
heat  transfer.  Furthermore,  a  large  specific  heat  per  unit  mass 
further  reduces  the  heat  transfer.  In  addition,  it  was  shown  that 
whereas  it  had  been  indicated  that  the  smaller  the  molecular  weight 
of  a  gas,  the  greater  its  effectiveness  in  an  injection  system, 
there  exists  the  possibility  that  the  combination  of  large  molecular 
diameter  and  a  large  molar  specific  heat  associated  with  polyatomic 
gases  may  in  fact  be  more  significant  in  reducing  skin  friction  and 
heat  transfer  than  the  low  molecular  weight  of  a  monatomic  gas  like 
helium.  (The  criterion  for  the  effectiveness  of  a  given  gas  is 
normally  established  as  the  relative  reduction  in  heat  transfer  for 
a  given  mass- injection  rate  of  the  gas.) 

In  the  present  investigation,  a  numerical  program  for  the 
solution  of  the  governing  equations  has  been  used  to  examine  the 
effects  of  helium,  air,  and  carbon-tetrachloride  injection  on  the 
skin  friction  and  heat  transfer  characteristics  of  the  lasinar 
boundary  layer  flow  over  a  flat  plate.  The  program  includes  the 
effects  of  equilibrium  vibrational  excitation,  that  is,  the  specific 


1 


NOLTR  66-215 


heat  is  assumed  to  vary  with  temperature .  Results  have  been  obtained 
for  Mach  numbers  from  0  to  3;  wall-to-stream  temperature  ratios  of 
0.25,  0.50,  and  1.00;  and  injection  rates  of  0,  -0.1,  and  -0.2.  The 
ranges  for  the  temperature  ratio  and  Mach  number  are  realistic  in 
terms  of  a  flat  plate  at  Mach  numbers  up  to  3,  and  for  the  forward 
portion  of  a  blunted  re-entry  body.  As  far  as  the  injection  rate 
is  concerned,  for  a  flat  plate,  separation  begins  at  a  value  of 
this  parameter  near  -0.6  according  to  our  definition.  For  very 
slender  bodies  (including  transverse  curvature  effects)  this  value 
may  be  as  high  as  -0.5  before  separation  takes  place.  For  an 
ablation  problem,  for  example,  a  slightly  blunted  cone  flying  at 
Mach  15  at  100,000  feet,  one  might  expect  values  of  the  injection 
rate  on  the  order  of  -0.4  to  -0.6.  Dissociation,  ionization  and 
thermal  diffusion  effects  have  been  neglected.  (A  study  of  the 
latter  appears  in  reference  (5).) 


ANALYSIS 


Basic  Equations  and  Boundary  Conditions 


We  consider  the  steady  two-dimensional  flow  of  a  pure  gas  forced 
into  the  stream  through  pores  of  the  flat  plate.  The  boundary  layer 
equations  describing  the  flow  of  the  binary  mixture  are  as  follows: 


Continuity: 


<pu>  +  -§^  (pv)  -  0 


(1) 


Momentum  in 
the  x 
direction: 


__  du 
9  "5y 


d  (a 

"3y  3y/ 


(2) 


Momentum  in 
the  y 
direction: 


(3) 


Conservation 
of  matter  of 
foreign  gas: 


xjW  of 

p*  lx  +  pv  -Sy 


+  aw(l-w) 


Energy: 


3b 


puT* 


State : 


p  —  pRT  **  constant 


(6) 
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In  these  equations,  u,  v,  p,  T,  h,  p,  and  k  represent  the  velocity 
coapcnents  along  and  noraal  to  the  plate,  density,  teaperature, 
enthalpy,  viscosity  coefficient,  and  coefficient  of  theraal  con¬ 
ductivity,  respectively.  In  addition,  w  is  the  mass  concentration 
of  the  foreign  gas,  defined  by: 


*lnl 


a,n,  -r  ^n. 


al“l 


(7) 


where : 


”  aoiecular  weight  of  species  "i" 
ni  ”  nuaber  of  aolecules  of  species  "i”  per  unit  voluae 
i  —  1  for  the  injected  gas 
i  *2  for  the  aain  streaa  gas 

Moreover,  in  equations  (l)-(6),  the  binary  diffusion  coefficient  is 
given  by  i,  the  thersal  diffusion  coefficient  by  k_,  and  the  theraal 
factor  by  a,  where:  1 


-  ow(l-w)  (S) 

The  gas  constant  H  is  for  the  aixture  per  unit  Bass. 

Equations  (1)— (6)  are  subject  to  the  following  boundary  conditions 


u<x,0)  -  G  (9) 

u(x,»)  -  u„  (10) 

T(x,»)  -  T.  (11) 

w(x,->  -  0  (12) 


v(2,0)  and  T(x,0)  prescribed 

It  should  be  pointed  out  that  the  iesperature  of  the  injected 
fluid  is  assused  identical  to  the  tesperature  of  the  plate  surface. 
No  consideration  is  given  to  the  source  of  fluid  nor  how  it  attains 
plate  surface  iesperature . 


Transformations 


As  Shea  (ref.  (3))  has  shown,  the  above  equations  and  boundary 
conditions  will  lead  to  sisilar  solutions  for  an  injection  rate 
varying  as  l//x  when  the  following  transformations  are  lutroduced: 
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and  the  stream  function: 

Y  -  (posUo^looX)  *  f  (t|) 


(13) 


(14) 


Under  these  transformations  then: 


pu 

ay 

1 

m 

~5y  s  p 

*5ti 

(15) 

pv 

-3Y 

1 

r  p»u 

af 

"  2 

3q  “  f) 

(16) 

pu 

d 

*55E 

+  pv 

5 

Hy 

i  -i 

-  PooUcoX 

f  d 

f  Sy  i 

(17) 

and  for  any  quantity  A(q)  we  have: 

~3y  [A(T,)ly]  “  1  p*u«x  1  TTn  {A(T1>  Hq}  (18> 


Next  we  introduce  nondiraensionaliaed  variables  (indicated  by 
barred  symbols)  by  dividing  all  flow  variables  by  their  corresponding 
stream  value,  except  for  w,  which  is  obtained  by  dividing  by  the  wall 
value  ww.  (Note  that  w  in  itself  is  already  dimensionless.  The 
introduction  of  w  is  for  the  convenience  of  certain  special  cases.) 


Using  then  equations  (15) -(18),  in  terms  of  the  dimensionless 
variables, equations  (l)-(6)  are  reduced  to  the  following  ordinary 
equations  obtained  by  Shen: 


d 

f—  d2f1 

!  f  d2f 

.  A 

HT 

I  +  f  — 2 

dr 

■  0 

d 

H? 

fag} 

,  f  dw 

+  f  w 

-1?  !>r] 

d 

nr 

[»  #] 

_  ds 

+  ' 

■I?  [?,iG] 

(19) 

(20) 

(21) 
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where : 

§ 

F 

G 


rn  -  rn  dTi 

new  independent  variable  -  j  pdq  -  J  ^ 

Cl  -  l]  f  +  a  I  (1-w*5)  h  lnT 


1  +  -££ 


fl_  p  dT  dE 
IPr  “  d? 


r  5  [' — fe-  —If  +  It  lnT]} 


(22) 

(23) 

(24) 


s  - 

h  -  \  u2  "  staSnation  enthalpy 

(25) 

Pi  - 

PC 

*  Prandtl  number  of  local  mixture 

uS  “ 

•  Schmidt  number  of  local  mixture 

Y  “ 

CP/CV 

In  deriving  equation  (24) ,  Shen  assumed  Cp  to  be  concentration 
dependent  [Cp(w)]  but  independent  of  temperature.  In  order  to 
include  temperature  dependence  in  Cp,  the  following  equation  for  G 
is  used  in  the  present  report  instead  of  (24) : 


G 


h2) 


(26) 


where  Le  -  Pr/w  «  Lewis  number.  This  formula  reduces  to  (24)  when 
a  =  0  and  Cp  is  assumed  independent  of  temperature. 

Now,  under  the  previous  transformations,  the  boundary  conditions 
(9)- (12)  become: 

(!)  |||  -  2  (ii)  |||  -  0 

U*>  ^0)  V*>  „Q 

(27) 

(iii)  f(0)  -  -2pv  /Re~ 

?=o  x 
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(i) 

w(0) 

-  1 

(ii) 

w  ( ®)  -  0 

(28) 

(i) 

S(0) 

*  Sw  -  constant 

(ii) 

S(«)  -  1 

(29) 

where  Rex  is  the  Reynolds  number  based  on  x  and  the  main  stream 
conditions,  “  ueox/v®. 

The  equations  describing  the  problem,  solved  numerically  in  the 
present  investigation,  are  equations  (19)- (21),  with  boundary 
conditions  (27)-(29). 


3.  Expressions  for  Transport  Properties 


A.  Binary  Diffusion  Coefficient 


From  reference  (6) ,  the  expression  for  the  binary  diffusion 
coefficient  in  practical  units  is: 


262.3  x 


i0“5 


[T3^ 


m.)/2  m.m  ] 
_ J. _ _ 


Pd^  nW 


(30) 


where : 


2 

diffusion  coefficient.,  cm  /sec 

•i(d.  +  d  .)>  angstroms,  (the  symbol  d.  stands  for 
*  ^  j  — 

molecular  diameter  or  collision  cross-section 
of  species  i) 


p  •  pressure,  atm. 


n 


T 

(1,D* 


temperature,  °K 

"reduced"  collision  integral  -  1  for  rigid  sphere 
model  as  used  in  the  present  investigation 


B.  Coefficient  of  Viscosity 

For  a  pure  gas  of  species  i,  reference  (6)  gives: 


^i 


266.93  x  10 


-7 


/m-jT 

12727* 


n' 


(31) 
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C.  Thermal  Conductivity 

For  a  monatomic  gas  of  species  i,  the  Chapman-Enskog 
Theory  yields: 


where : 


k .  - 


R,  - 


k,  «  ~  R.g. 
i  4  11 

thermal  conductivity,  cal/cra-sec-°K 
viscosity  given  by  equation  (31) 
gas  constant  for  species  i 


(37) 


Substitution  of  (31)  into  (37)  yields: 


-4  1 

k_.  -  1.9889  x  10  / Y 


Nondimensionalizing  (38)  we  have: 

-  kl 

kl  “  kt  “ 
k2 

k2  -  ir  “ 


mi  di2 


#$)  * 


(33) 


(39) 


(40) 


For  a  mixture  of  monatomic  gases,  reference  (6)  has: 


1  +  (1.065)  G12^2  1  +  (1.065)  G21fl 

nl  n2 


(41) 


The  above  formulas  can  be  extended  to  polyatomic  gases,  as 
indicated  by  Dorrance  in  reference  (6) .  The  correction  factor  Eu 
(the  Eucken  correction  factor)  is: 


Eu  - 
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where : 

k'  -  corrected  thermal  conductivity  for  polyatomic  gases 

k  -  thermal  conductivity  calculated  from  monatomic 
considerations 

,  pS. (5/2)1^ 

L  -  - £ - 

R  -  R^iik  =»  universal  gas  constant 

L'  can  be  evaluated  using  the  Chapman-Enskog  theory  for 
monatomic  molecules.  As  indicated  in  reference  (6),  a  good  value 
is  L'  s  0.885.  Using  this  value,  the  correction  factor  becomes: 

Eu  -  kVk  -  0.115  +  0.354  (42) 


Furthermore,  the  Chapman-Enskog  expansion  for  the  thermal  conductivity 
of  a  mixture  of  monatomic  gases  can  be  modified  making  use  of  the 
Eucken  factor  to  yield  an  expression  for  the  thermal  conductivity  of 
polyatomic  gases.  The  equation  is  (for  a  binary  mixture) : 


1  +  1.065  G12“2 
nl 


kl  ^U1 


1  +  1.065  G12_2 
nl 


1  +  1.065  G, 


21— 


k2  Eu2 


1  +  1.065  G21^l 
n2 


Nondimensionalizing,  we  obtain: 


k 


-  k/k 

CD 


5,  Eu, 

1  X _ 

1  +  1.065  G 


n2 

12-~ 

nl 


+ 


g2  £u2 _ 

1  +  1.065  G21jL 
n2 


(43) 


where  and  £3  are  given  by  (39)  and  (40),  respectively,  and  Eu  by 
(42).  In  addition,  G^  is  given  by  (33)  for  -  T,,. 
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4.  Enthalpy  and  Specific  Heat 

The  total  internal  energy,  E,  of  one  mole  of  a  perfect  gas, 
including  translational  as  well  as  internal  degrees  of  freedom, 
may  be  expressed  by  the  following  equation  as  indicated  by  Herzberg 
in  reference  (7) : 


E 


i 


V2  HT  <Xn<ii* 


(44) 


where  the  subscript  i  indicates  that  we  are  talking  about  a  given 
gas  of  ith  species.  In  this  equation,  Q  is  the  total  partition 
function,  E°  the  energy  at  absolute  zero  temperature,  T  the  absolute 
temperature,  and  K  the  gas  constant. 

Now,  the  enthalpy,  H,  of  one  mole  of  a  perfect  gas  is  the  sum 
of  the  total  internal  energy,  E,  and  the  external  energy  pv  *  RT, 
such  that: 


H . 


-  +  R±T  + 


RiT2  <t«Qt) 


(45) 


Finally,  the  specific  heat  at  constant  pressure  is  given  by: 


dHi 

cIT" 


Ri  +  Ri  HT 


[t2  It  (W 


(46) 


Now,  Q  can  always  be  written  as  the  product  of  the  translational 
and  internal  partition  function  QxRi  and  ^INT*  Therefore,  since  Q 
occurs  in  (45)  and  (46)  only  as  fnQ,  then  the  enthalpy  and  specific  heat 
are  the  suras  of  a  "translational”  and  an  "internal"  term.  In 
particular,  for: 


S'R 


1.879  x  1020 


3/2 


V  IB 


(where  V  is  the  volume  considered,  ra  the  chemical  atomic  or  molecular 
weight,  T  the  temperature. . .for  details  see  Herzberg,  reference  (7)), 
we  can  write: 


H 


i 


5 

1 


+  1  RiT  +  HINT 

(47) 

+  Cp 

(48) 

HINTi 
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where : 


HINTi 

cn 

PlN’T 


RiT2  ^  (*hQTMfr) 


Ei?r  [• 


INT 

2  d 

3T  j 


(49) 


Moreover,  can  be  v/ritten  as  the  product  of  a  number  of 

factors  if  certain  approximations  are  made.  In  particular,  if  the 
interaction  of  vibration  and  rotation  is  ueglected,  both  hIN  -j-  and 
Cpi^r  can  be  written  as  the  sums  of  a  rotational  and  vibrational 

term  such  that: 


liINTjl 


where : 

HROT . 

1 

C 

pROTi 

HVIB. 

c 

?VIB . 

l 


iIROTi  +  HVIBi 

»  c  +  c 

pROT .  PVIB . 

l  x 

Ri  T  ~dT  (XnQROT) 

Ri  "ar  [T  dT  (XiiQrot] 
Ri  T  cTT  ^in<*VIB^ 

"  Ri  It  [t2  It  (£nQviB}] 


A.  Monatomic  Molecules:  Helium 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


For  monatomic  molecules,  in  which  only  the  translational 
degrees  of  freedom  exist,  we  have  from  (43) : 


5  R 
2  Ri 


(56) 
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and  from  (47) ,  for  E? 


0  we  have : 


IF 


I  Ri  T 


Cp.  T 
*1 


(57) 


For  helium  R-  -  0.4967  cal/gm-  K  such  that  Cp.  »  1.24 

o  XHe 

cal/gm-  K.  These  are  the  values  used  in  the  present  investigation. 


B.  Diatomic  Molecules 


(1)  Vibrational  Excitation  Neglected 

The  effects  of  translation  offer  nothing  new,  being 
the  same  for  diatomic  as  well  as  for  monatomic  gases  (as  well  as 
for  polyatomic) . 

The  equations  for  Qrot  for  molecules  without  free  or 
hindered  internal  rotations  are  given  by  Herzberg  in  reference  (7) . 
From  the  equation  for  the  case  of  diatomic  (and  linea/i.-  polyatomic) 
molecules  we  have,  after  substituting  in  (52)  and  (53): 


H 


ROT. 


Ri  T 


Therefore,  neglecting  the  vibrational  contribution  we 
have,  for  diatomic  molecules: 


H.  -  I  R ,  T  -  C  T 
i  2  i 

cPi  ' 


For  air,  then,  considered  as  a  "mythical  diatomic 
molecule"  we  have  R.,  —  0.0685  cal/gm-°K  and  thus  C„  -  0.24 

iAIR 

cal/gm-°K.  These  are  the  values  used  in  the  present  investigation. 


(2)  Vibrational  Excitation  Included 

For  diatomic  molecules  with  vibrational  excitation 
included,  the  Hjifr  term  in  equation  (47)  contains  the  contribution 
due  to  rotation,  T  as  before,  plus  now  a  vibrational  term  coming 
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from  Qvib  which  we  shall  denote  as  Ev.  This  "vibrational  energy"  can 
be  written  as  (see  reference  (8)): 


R 


i 


-1 


(60) 


where : 

9  «  Planck’s  constant  x  vibration  frequency 

v  Boltzmann's  constant 

•  characteristic  vibrational  temperature  for  given 
gas 


vibrational  temperarure  (for  equilibrium  -  T) 


Hence,  the  expression  for  the  enthalpy  is,  assuming  equilibrium: 


Hi  "  |  Ri  T  +  Ey  (T)  (61) 

and  the  specific  heat  at  constant  pressure  is: 

da<  7  7 

CP.  <T>  -  <rr--2Ri +  CP.  m  <62) 

^VIB 

where  the  prime  denotes  differentiation  with  respect  to  T  .  In 
particular ,  from  equation  (60) , 

By  /2T 

V  m  -  Ri  {— h-8  -■«}  (S3) 

sxnh  vi/ 

In  the  present  investigation,  as  in  reference  (9),  air  is  treated  as 
a  "mythical  diatomic"  gas  with  Ri  -  0.06851  cal/gm-°K  and  8y  =*  3076°K. 


C.  Polyatomic  Molecules 


In  this  case,  the  effects  of  translation  offer  nothing  new, 
being  the  same  as  for  the  monatomic  and  diatomic  gases.  For  linear 
molecules,  the  treatment  of  rotation  given  for  diatomic  gases  is 
applicable.  For  nonlinear  molecules,  such  as  CCL4  and  CH4,  the 
effect  of  rotation  adds  (3/2)R^T  to  the  enthalpy,  such  that: 


Hi _  -  %  R4T  -  Cn  T 

2  i  ProT 


lROT 


(64) 
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There  are,  on  the  other  hand,  several  modes  of  vibration, 
some  of  these  important  at  low  temperatures,  others  at  high  tempera¬ 
tures.  In  addition,  there  are  degeneracies  in  these  nodes  so  that 
one  has  less  "fundamental"  modes  available. 

However,  using  the  "harmonic-oscillator"  approximation  (see 
reference  (7)),  Qyiy  is  a  product  of  terms  due  to  different  normal 
vibrations.  The  resulting  expression  for  llyjg  then  turns  out  to  be: 


II.  -  ^  E 

XVIB  tT 


(65) 


where  indicates  the  summation  over  all  fundamentals  of  the  molecule 
and: 


where : 


•« 


n 


Ri 


e 


3 — 7T 
v. 


1  -1 


2.858 

“TFT  '°n 

frequency  of  vibration  in  cm' 


(66) 


(67) 


d 


n 


degrees  of  degeneracy  in  vibration  corresponding 
to  >jb 

n 


That  is  to  say,  for  these  polyatomic  molecules,  the  vibra¬ 
tional  energy  term  Ev  (T)  encountered  in  the  discussion  of  diatonic 
molecules  now  consists  of  a  summation  of  "vibrational  energy  terms” 
corresponding  to  the  different  fundamental  frequencies  of  vibration, 
multiplied  by  the  degrees  of  degeneracy  in  each. 


For  example,  methane  has  four  fundamental  frequencies  of 
vibration,  i.e.,  n  *  1,  2,  3,  4.  Therefore: 


H 


VIB 


CH, 


E  (T)  -  E  -f  E  -r  E  +  E 

v  V1  v2  v3  1 


2.353 


R 


R„,t  {d 


CH, 


ch*  n _ j2 ,35ti  r  - 

4  exp^Tr  lJ 

4 


-}-  term  for  n~2 
+  term  for  n«3 
-f  term  for  n-4 
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In  particular,  from  reference  (7): 
0)1  “  2914.2  era"**  '«2 
d.  -  1  d2  * 


1526  ca 
2 


3032.3  cm 
3 

RCH  "  °*1242 

n4 


m  1306.2  ca 
d4  -  3 

cal/cm-°K 


(molecular  weight  of  CH^  —  16} 


The  specific 


c  (t; 

VIB 

which  is,  except  for  the  summation  sign  and  the  factor  dfi.  the  same 
as  (62)  and  (63)  for  diatomic  molecules.  That  is,  the  vibrational 
contribution  Co  the  specific  heat  can  also  be  represented  as  the  sua 
of  smaller  "specific  neats’*  corresponding  to  the  fundamental  modes 
of  vibration  multiplied  by  the  degrees  of  degeneracy  in  each. 


heat  at  constant  pressure  is  given  by: 


(T) 


Ri  L  dn 

a 


I  V» 

isiilh 


V/2T  -j2 

(9;  /T)i 

n 


(68) 


then: 


or 


Finally,  the  total  enthalpy  for  these  polyatomic  gases  is 


H.  -  5/2  R|T 

(translation) 


3/2  EjT 
(rotation) 


Ye 

4»  v 
n  n 

(vibration, 
given  by  (66) ) 


H. 


l 


4R±T 


(69) 


The  specific 


heat  at  constant  pressure  is: 
-  5/2  Rj  4-  3/2 

(translation)  (rotation) 


(vibration, 
given  by  (53)) 


(70) 
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If  we  denote  by  C  the  contributions  to  the  specific 

PT+U 

heat  vvliich  are  independent  of  temperature  and  thus  constant  for  a 
given  gas,  i.e.,  translation  and  rotation,  then  we  have  the  more 
convenient  expressions: 


II. 


C  .  T+  E,  (T) 
‘"f+U  V 


(71) 


C 

pT+it 


4-  li 


(T) 


(72) 


.vhich  apply  to  all  ga-ies  (with  the  values  of  C 
different  for  each,  of  course). 


and  E  (T)  being 
pT+il  V 


Suppose  all  the  vibrational  moaes  of  a  given  gas  are  excited, 
say  CCL^  (as  in  reference  (4)).  Then  for  each  vibrational  degree  of 

freedom  th *  contribution  to  C .  is  Hi .  Therefore,  fo~  the  4.7 

PVIE 

vibrational  degrees  of  freedom  in  CCi,^  O'-CCL^  =  .0129  cal/gm-°K) 

we  have:  C  -  4.7  (.0129)  *  .0604  cal/gm-°K.  The  total  specific 
PVIB 

heat  at  constant  pressure  is  then: 


411 


CCL , 


+  cp 


VIB 


(translation  +  rotation) 


.0516  -f  .0604 
.1120  cal/gm-°K 


The  specific  heat  at  constant  pressure  of  the  mixture  is 
either  concentra cion-dependent  alone  (vibrational  excitation 
neglected) ,  or  both  concentration  and  temperature-dependent  (vibra¬ 
tional  excitation  included).  In  any  case, 


(1*  w) 


(73) 


where  C  and  C,  are  either  C  or  0 

P1  p2  PT+U  PT+R 

whether  vibrational  excitation  is  neglected 


+  C  depending  on 
PVIB 

or  included > 
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Nondimensionalizing  (73),  we  have: 


w  C  +  (1-w)  C 


P« 


(74) 


where  in  the  present  investigation  we  have  considered 


(CP.  >2  Tc 

hT+R 


(75) 


always.  The  same  reasoning  applies  for  the  enthalpy  of  the  mixture: 


H 


w  H1  +  (1-w)  H2 


ar 


w  H1  +  (1-w)  H, 

3 T-  ‘ 

ft+r 


(76) 


where  and  H2  either  neglect  or  include  vibrational  excitation. 

5.  Energy  Transfer  Across  Interface 

As  indicated  by  Korobkin  in  reference  (4) ,  the  "energy  flux 
across  the  interface"  (i.e.,  the  energy  which  will  increase  the 
body  temperature)  can  be  expressed  as: 


<-E)w  -  [k  If  -  VT]W 


(78) 


Usin[  the  transformations  outlined  previously,  this  can  be 
written  iu  the  dimensionless  form: 


(-E)w  -  Ew  Pw 


-  T 


w 


w 


f(0) 


Pr 


CPl 

CP2 


(79) 


That  is,  the  surface  heat  transfer  is  seen  to  be  divided  into 
contributions  coming  from  a  "temperature  gradient"  term  and  a 
"concentration  term.' 


by: 


Now,  it  is  desired  to  reduce  (-E)^.  Evidently,  this  can  be  done 


a.  Decreasing  k  .  This  may  be  accomplished,  for  example,  by 
injecting  a  large  molecular  diameter  gas,  since  £  — ~  . 

w  d| 
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b.  Increasing  T  .  This  will  increase  the  magnitude  of  the 

w 

negative  term. 

(C  /mole) 

c  Increasing  Cpj.  Since  -  ■^olVcuTar  wglgKtTy  > 

this  may  be  accomplished  by  either  having  an  injectant  of  small 
molecular  weight,  or  one  with  a  high  specific  heat  per  mole  (say, 
as  a  result  of  vibrational  excitation) . 

Evidently,  the  relative  efficiency  of  any  of  these  actions 
depends  on  the  importance  of  the  term  in  which  the  affected 
parameter  appears.  For  example,  if  the  injection  is  going  to  be 
small  (small  f(0))  then  perhaps  the  heat  reduction  efforts  should 
be  concentrated  on  the  first  terra  of  (79),  i.e.,  on  reducing  Ew, 
whereas  if  the  injection  is  going  to  be  considerable,  the  second 
term  might  be  more  important. 

In  addition,  all  of  the  factors,  (a),  (b) ,  (c) ,  do  not  necessarily 
go  together.  And  so  one  might  have  a  large  molecular  diameter  gas 
but  with  a  high  molecular  weight  (an  effect  which  may  in  turn  be 
corrected  by  a  high  specific  heat  per  unit  mass),  and  so  on. 

In  sum,  the  ideal  situation  is  to  have  a  favorable  combination 
of  all  these  factors.  Reflecting  this  fact,  we  will  see  in  the 
present  report  that  at  some  wall-to-free  stream  temperature  ratios, 
the  light  gas  (helium  in  our  case)  is  the  most  effective  coolant, 
whereas  at  other  temperature  ratios,  the  heavier  gas  is  more  advan¬ 
tageous. 


6.  Numerical  Solution  of  the  Equations 


Equations  (19)-(21),  subject  to  the  boundary  conditions  (27)-(29), 
have  been  solved  on  an  IBM  7090  computer  using  the  program  known  as 
FNOL  2  described  in  reference  (10).  This  program  is  a  subroutine 
which  employs  the  Runge-Kutta  and  Adams-Mouiton  fourth  order  tech¬ 
niques  for  the  integration  of  a  system  of  ordinary  differential 
equations  (up  to  30). 

Mathematically,  the  problem  consists  in  solving  three  ordinary 
differential  equations  for  the  three  functions  f(0,  s(§),  and 
w(s).  These  equations  are  linked.  In  order  to  use  FNOL  2  ,  these 
three  equations  are  reduced  to  a  system  of  seven  first  order  ordinary 
differential  equations,  given  by: 

D.  -  dY*  -  zH  v  where  Y.  -  pp  (30a) 

1  <*§  Y1  1  d52 


18 


NOLTR  66-215 


Equations  (80a) , 


dY2 

nr 

m 

•i-  Y.  where 
PP 

Y  - 

x2  d5 

(80b) 

dY3 

nr 

- 

Yg  where 

Y3  "  f 

(80c) 

dY4 

w~ 

-  - 

Y.  where 
W  4 

Y4  - 

u> 

(30d) 

dYS 

d? 

- 

m 

—  Y4  where 
PP 

y5  -  i 

(80e) 

dY6 

nr 

-  - 

f  ^  where 

Y6  ”  w  (§f +  G) 

!  (30f ) 

dY? 

nr 

- 

^  where 

Y?  -  s 

(80g) 

(80b) 

,  and 

(80c)  come  from  (19);  equations 

(80d) , 

and  (80e)  come  from  (20) ;  and  equations  (30f )  and  (80g)  from  (21) . 


Now,  the  following  boundary  conditions  must  be  given  for  solution 
of  these  equations:  f(0);  f'(0),  f#(0),  s(0),  S'(0),  w(0),  w'(0) 
where  the  prime  denotes  differentiation  with  respect  to  5.  Unfor¬ 
tunately,  not  all  of  these  conditions  are  actually  given.  In 
particular,  we  are  missing  f'(0),  s#(0)  and  w'(0).  But  we  do  have 
f#(®),  s(«*),  and  w(®).  Therefore,  a  scheme  must  be  devised  to 
obtain  the  missing  boundary  condi tiors  which  will  result  in  given 
conditions  at  §  *  •.  Once  they  arc  .btained,  the  differential 
equations  may  be  solved  numerically  using  FN0L2. 


The  procedure  to  obtain  the  required  boundary  conditions  used  in 
the  present  investigation  is  based  on  the  one  employed  in  ref¬ 
erence  (11).  Stated  abstractly,  we  have  three  functions  (f,g,h)  of 
three  variables  (x,y,z)  and  we  wish  to  find  the  point  at  which  they 
are  zero.  In  particular,  if 

f (x,y,z)  -  f'(«*>  -  2 

g(x,y,z)  -  s(»)  -  1 

h(x,y,z)  -  w(«) 


where  x  *  f'(0),  y  -  I'(0)  and  z  -  w(0),  _we  ask  for  what  values  of 
x,  y  and  z  are  f-g-h-0  (i.e.,  f'(«)  «  2,  s'(*)  -  1,  w  (•)  «  o) . 
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Such  a  problem  can  be  solved  using  the  Newton-Raphson  iteration 
method.  The  iteration  equations  are,  in  matrix  notation: 


W  tl 


rxx 

(yj 


z  n 


df  at  af 

"Sx  "Sy  ~Sz 

dg  dg  ag 

"ox  "cfy  "5z 

Sh  dh  dh 
\  "5x  "Sy  "5z 


(81) 


This  iteration  would  then  be  accomplished  as  follows: 

(a)  Assume  first  trial  values  of  f'(0),  s#(0)  and  w(0).  (These 
are  x0,  yQ,  z0,  respectively.)  Integrate  the  differential  equations 
(30)  using  FNOL2  obtaining  values  of  f#(«),  s(«*),  w(®).  These 
are  fQ,  go»  ho*  respectively. 

(b)  Three  more  complete  integrations  are  made,  each  differing 
from  that  of  (a)  bg  a  small  perturbation  of  one  of  the  initial 
conditions  f'^O),  s'(0),  w(0),  i.e.,  (x,y,z).  In  each  case,  values 
of  f/(®),  s(«>)  and  w(«»)  are  obtained. 

(c)  From  the  results  of  (a)  and  (b)  approximate  values  of 

etc.  are  computed  in  an  obvious  manner. 

(d)  The  matrix  of  partial  derivatives  is  inverted  and  new 
values  of  x,  y,  z  are  computed  from  (31). 

(e)  The  process  stops  when  f,  g,  h  are  as  close  to  zero  as 
specified  (depending  on  the  accuracy  desired).  The  last  x,  y,  z 
are  the  missing  values  of  f * (0) ,  s#(0)  and  w(0)  which  meet  the 
known  boundary  conditions  at  infinity,  and  the  Integration  with 
these  values  is  the  solution  to  (80)  with  the  appropriate  boundary 
conditions. 

Additional  details  are  given  in  reference  (y) .  The  interval  of 
integration  used  in  the  present  investigation  was  As  =*  0.1  for  the 
case  of  no  vibrational  excitation,  and  As  =*  0.01  with  vibrational 
effects  included. 


RESULTS  AND  DISCUSSION 

As  explained  in  the  Introduction,  the  main  purpose  of  this 
part  of  our  investigation  was  to  obtain  boundary  layer  profiles 
for  evaluation  of  their  stability  characteristics  by  the  i_ethoa 
described  in  reference  (1) .  In  addition,  once  the  profiles  were 
obtained,  the  effect  of  foreign  gas  injection  on  their  skin  friction 
and  heat  transfer  characteristics  was  examined. 
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Figure  (1)  presents  the  effects  of  injection  on  the  skin  friction 
coefficient  at  a  representative  Mach  number  (0.70)  and  wall-to- 
stream  temperature  ratio  (0.50).  The  expression  for  C^/llex  is 
obtained  from: 


C 


p®[du/dyj0 

1  ^ 

■k  P  a 

OO  OD 


which  becomes,  after  the  transformations  have  been  applied: 


C£^ 


2^-s 


fl  —  d~fl 
7  PP  — n 


(32) 


We  see  in  the  figure  that  helium,  the  lightest  gas,  gives  the 
sharpest  reduction  in  the  skin  friction  coefficient  as  the  injection 
is  increased.  This  is  due  to_the  effect  of  molecular  weight  of  the 
injected  gas  on  the  quantity  pjl.  In  figure  (2)  the  injection  rate 
is  held  constant  at  the  same  Mach  number,  and  the  wall-to-stream 
temperature  ratio  is  varied.  Again  helium  appears  most  favorable. 
Finally,  in  figure  (3)  the  Mach  number  is  varied  and  the  injection 
rate  and  temperature  ratio  are  held  constant  with  similar  results. 

So  we  see  that  for  these  ranges  of  Mach  number,  temperature  ratio, 
and  injection  rate  the  effect  of  injection  on  the  skin  friction  is 
directly  related  to  molecular  weight. 

The  heat  transfer  results  are  not  as  straightforward  as  in  the 
case  of  skin  friction.  Consider  figure  (4)  which  shows  the  effects 
of  wall-to-stream  temperature  ratio  on  the  ratio  of  Stanton  number 
with  injection  to  Stanton  number  with  no  injection.  The  figure 
might  best  be  understood  by  looking  at  equation  (79) .  We  see  that 
at  temperature  ratios  greater  than  .36,  helium,  the  lightest  gas, 
appears  as  the  most  effective  coolant.  This  is  because  the  value 
of  Tw  is  high  enough  to  make  important  the  second  term  in  (79) . 
Therefore,  the  reduction  in  heat  transfer  achieved  by  injecting  a 
small  molecular  weight  gas  (which  further  increases  this  term)  is 
further  enhanced  by  the  large  value  of  Tw.  However,  at  lower  values 
of  Tw,  this  term  is  no  longer  important,  aDd  effort  should  be  con¬ 
centrated  in  decreasing  the  first  terra.  This  may  be  accomplished 
by  injecting  a  gas  with  large  molecular  diameter,  and  thus  decreasing 
icw.  The  figure  indeed  shows  that  at  lower  temperature  ratios,  the 
heavier  gases  are  more  effective  than  helium  in  reducing  heat 
transfer. 

In  figure  (5) ,  the  results  are  plotted  versus  Mach  number  for 
various  temperature  ratios  and  an  injection  rate  of -9.1.  In  addition 
to  confirming  again  the  preceding  analysis,  the  figure  shows  that 
the  results  for  a  given  Tw  do  not  vary  much  with  Mach  number  (except 
Tw  -  0.50  for  helium  injection).  In  fact,  a  correlation  has  been 
suggested  by  Gross  and  others  (ref.  (12))  to  collapse  all  the  data  at 
different  Mach  numbers  and  temperature  ratios  into  one  curve.  This  is 
done  by  plotting  the  Stanton  number  ratio  versus  one  half  the  injection 
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FIG.  1  EFFECT  OF  INJECTION  RATE  ON  SKIN  FRICTION  WITH 
LIGHT,  MEDIUM,  HEAVY  INJECTANTS 
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FIG.  2  EFFECT  OF  WAIL  -  STREAM  TEMPERATURE  RATIO  ON  SKIN 
FRICTION  WITH  LIGHT,  MEDIUM,  HEAVY  INJECTANTS 
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I  {  ( O }  =  -0.10 

j  Tw/T-  =  0.50 

0.1 1- 

! 

S 


MACH  NUMBER 


FIG.  3  MACH  NUMBER  EFFECT  ON  SKIN  FRICTION  WITH  LIGHT, 
MEDIUM,  HEAVY  INJECTANTS 
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rate  parameter  divided  by  the  square  root  of  the  Chapaan-Kubesin 
constant  evaluated  at  a  reference  temperature  (C*) .  This  is  pre¬ 
sented  in  figure  (6) .  The  figure  shows  that  the  correlation  worked 
for  the  Mach  number  range  considered  (Mach  numbers  of  0  to  3)  but 
failed  as  far  as  the  temperature  ratio  was  concerned.  Again  the 
results  predicted  by  equation  (79)  appear:  whereas  at  a  temperature 
ratio  of  0.5,  helum  is  the  most  effective  coolant;  when  the  tempera¬ 
ture  ratio  is  as  low  as  0.25,  the  heavier  gases  are  more  advantageous. 

We  have  seen  then  that  there  are  conditions  when  it  is  more 
favorable  to  use  a  heavy  injected  gas  instead  of  having  a  light  gas 
injection  system.  This  is  even  more  so  when  one  considers  the 
effect  of  injection  cm  the  stability  characteristics  of  the  boundary 
layer  as  in  reference  (1) .  Here  it  is  shown  that  the  heavy  gas 
can  actually  cause  the  boundary  layer  to  be  more  stable,  whereas 
a  light  gas  is  destabilizing. 


CONCLUSION 

The  main  purpose  of  the  investigation  was  to  obtain  boundary 
layer  profiles  showing  the  effects  of  a  light,  medium,  and  heavy 
gas  injection.  The  method  of  solution  of  tLe  governing  equations 
has  been  described.  The  stability  characteristics  of  the  resulting 
profiles  appear  in  reference  (1).  In  addition,  the  effects  of 
injection  on  the  skin  friction  and  heat  transfer  characteristics 
of  these  profiles  have  been  obtained.  The  results  have  shown  that 
the  reduction  of  skin  frictiou  and  heat  transfer  depends  in  general 
on  factors  which  are  not  necessarily  complimentary .  For  exaaple, 
sometimes  a  large  molecular  diameter  is  advantageous,  other  times  a 
small  molecular  weight  is  better,  and  so  on.  Furthermore,  we  have 
seen  that  the  value  of  the  wall-to-streae  temperature  ratio  will 
also  affect  the  choice  of  Injectant .  In  short,  the  selection  of 
discrete  injection  systems  must  represent  a  favorable  compromise. 
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